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Abstract 



In this paper, the cosmic microwave background (CMB) anisotropy in a 
multiply-connected compact flat 3-torus model with the cosmological constant 
is investigated. Using the COBE-DMR 4- year data, a full Bayesian analysis 
revealed that the constraint on the topology of the flat 3-torus model with 
low-matter-density is less stringent. As in compact hyperbolic models, the 
large-angle temperature fluctuations can be produced as the gravitational 
potential decays at the A-dominant epoch well after the last scattering. The 
maximum allowed number N of images of the cell (fundamental domain) 
within the observable region at present is approximately 49 for f2 m = 0.1 and 

= 0.9 whereas N ~ 8 for Q m = 1.0 and £7a = 0. 
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I. INTRODUCTION 



For a long time, cosmologists have assumed the simply connectivity of the spatial hy- 
persurface of the universe. If it is the case, the topology of closed 3-spaces is limited to 
that of a 3-sphere if Poincare's conjecture is correct. However, if we assume that the spatial 
hypersurface is multiply connected, then the geometry of spatially finite models can be flat 
or hyperbolic as well. The metric describes only the local geometry. We should be able to 
observe the imprint of the "finiteness" of the spatial geometry if it is multiply connected 
on scales of the order of the particle horizon or less, in other words, if we live in a "small 



universe" . 



For flat multiply connected models without the cosmo logical constant, various constraints 
using the COBE DMR data have been obtained |T]-0. Assuming that the initial power 
spectrum is scale-invariant (n—1) the suppression of the temperature fluctuations on scales 
beyond the typical size of the cell L[| leads to a decrease in the large-angle power. For 
3-torus models without the cosmological constant in which the cell (fundamental domain) 
is a cube, the constraint is L > 0.8-R* where R* is the comoving radius of the last scattering 
surface (horizon radius) [|^@]- It should be emphasised that the constraint itself does not 
imply that the "small universe" is ruled out since the possible maximum expected number 
of the copies of the cell within the observable region at present is approximately 8. For 
models in which one side of the cell is longer than the others, the constraint for the smallest 
topological identification scale(=the diameter of the smallest ball which can wrap around 
the space) can be less stringent 0. 

On the other hand, recent observations of distant supernova la imply the existence 
of "missing energy " which possesses a negative pressure and equation-of-state (w = p/ p) in 
the form of either a cosmological constant, vacuum energy density or a slowly- varying spa- 
tially inhomogeneous component "quintessence". Together with the first observation of the 



height and the position of the first Doppler peak by Boomerang |I| and MAXIMA |TT 



and the COBE-DMR data, the constraint for the cosmological constant is 0.69 < < 0.82 
and for the matter(cold dark matter, baryon plus radiation) 0.28 < Q m < 0.42 assuming 
adiabatic initial perturbations [IT2 . 



For low-matter-density models with flat geometry a bulk of large-angle cosmic microwave 
background (CMB) fluctuations can be produced as the gravitational potential decays at the 
A-dominant epoch 1 + z ~ Recent works |13|HT6"1| have shown that the angular 

power spectrum C\ is completely consistent with the COBE-DMR data for some compact 
hyperbolic models which are incompatible with the previous analyses JD]]. Because the an- 
gular sizes of fluctuations produced at the late epoch are large compared to those on the 
last scattering for flat or hyperbolic geometry, we expect that the constraints for compact 
flat models with low-matter-density can be also significantly loosened. 

In this paper, the CMB anisotropy in a multiply-connected compact flat 3-torus model 
with or without the cosmological constant is investigated. In sec II we briefly describe the 
time evolution of the scalar perturbation in the locally flat Friedmann-Robertson- Walker 



L can be defined as twice the diameter which is defined as the maximum of the minimum geodesic 
distance between two points over the space. 



2 



models which will be used for computing the CMB anisotropy. In sec III we compare the 
angular power spectrum in a flat 3-torus model with the cosmological constant to the one 
in the "standard" 3-torus model with Q tot = 1. In sec IV a full Bayesian analysis using the 
COBE-DMR data has been carried out for giving the constraint on the minimum size of the 
cell. 



II. CMB ANISOTROPY 

In what follows we assume that the matter consists of two components:a relativistic one 
and a non-relativistic one (fl m = Q r + First of all we consider the evolution of the 

background space. From the Friedmann equation, the integral representation for the time 
evolution of the scale factor a (normalised to 1 at present time) in terms of the conformal 
time 7] takes the form 

rj( a ) = H Q ~ l f da == (1) 

Jo vf2 n0 ci + ^rO + ^acl 

where H Q , Q rQ and f2 n0 are the Hubble parameter, the density parameter for relativistic and 
non-relativistic matter at present, respectively. Note that (|1|) can be written in terms of the 
elliptic integral of the first kind which is too complex to describe here. The radius of the 
last scattering surface R* in comoving coordinates can be obtained by integrating rj from 
1/( z lss) ~ 1/1100 to 1. One can easily see that the presence of Q\ makes R* to a larger 
value. For a fixed H and the total density at present VL tot) the presence of a cosmological 
constant or vacuum energy implies a decrease in the cosmological expansion rate in earlier 
epoch which augments the age of the universe and the size of the causally connected region. 
As for the radiation density, the standard scenario of the thermal history gives the present 
density Q r o = Q(photon) + Q(neutrino) = 4.15 x 10~ 5 /i~ 2 assuming three generations for 
neutrinos [18|]- Applying the value to ([!]) one obtains i?* = 4.98 for (Q n0 ,Q/C) = (0.1,0.9) 
while = 1.99 for (tt n0 ,tt A ) = (1.0,0). 

Next, we consider the time evolution of the scalar- type perturbation. Assuming that 
the anisotropic pressure is negligible, in the Newtonian gauge, the Newtonian curvature 



perturbation $ in k space satisfies [fLS 



+ 3H(1 + c 2 s )& + c 2 s k 2 $ + (2W + (1 + 3c 2 s )H 2 )<£ = —AnGa 2 pT, (2) 



2 _ 1 Pr 

3 Pn + \Pt 
p _ Pn ( 3 5pr_ _ Sj)n 

Pr + \pn V 4 Pr Pn 

where ' denotes the conformal time derivative d/drj, TC = a' /a, G is the Newton's constant, 
p is the total pressure, c s is the sound speed of the fluid and Y corresponds to the entropy 
perturbation defined as pY = 5p — ^<5p. Assuming a pure adiabatic perturbation (r = 0), 
the amplitudes of $ for the non-decaying mode are constant in the radiation dominant 
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epoch if krj eq « 1 Hence the non-decaying mode of the perturbation that enters the 
Hubble radius after the radiation-matter equality time t] eq can be obtained by setting the 
initial condition $'(0) = 0. In practice, we approximate the solution of ([I]) by interpolating 
polynomial functions which are used for numerically solving the second order ordinary dif- 
ferential equation (0). 

For superhorizon perturbations, the amplitudes of $ are constant in the radiation dom- 
inant epoch ($/$(0) = 1) and in the matter dominant epoch ($/$(0) = 9/10) However, in 
the A-dominant epoch $ gradually decays as 1/a that will play a critical role in character- 
ising the large-angle temperature fluctuations for A-models. 

Assuming a pure adiabatic initial condition and neglecting the anisotropic pressure, the 
temperature fluctuation on large angular scales can be written in terms of the curvature 
perturbation $ as 

AT. s 1 . . . rvo d®(rj, (ri — 77)11) , , . 

(n) = --$(//*, (770 - V» ~ 2 / 1 K '° ^77, (4) 



where n denotes the unit vector which points towards the sky and 77* and 770 correspond 
to the last scattering and the present conformal time, respectively |]20|j2T[| . The first term 
in the right-hand side in (^) describes the ordinary Sachs- Wolfe (OSW) effect while the 
second term describes the integrated Sachs- Wolfe (ISW) effect. Note that the Sachs- Wolfe 
formula (^) is valid only for perturbations whose scale is sufficiently larger than the horizon 
radius at the last scattering k < I/77*. From ([!]), one obtains the conditions k < 20.5H 
for (n m ,Q A ) = (1.0,0) and k < 9.9H for (n m ,Q A ) = (0.1,0.9). The effect of the acoustic 
oscillations becomes important for fluctuations on smaller scales. However, for large angular 
scales (/ < 15), contributions from acoustic oscillations are negligible. 



III. POWER SPECTRUM 

Let us now consider cosmological models whose spatial geometry is represented by a flat 
3-torus which is obtained by gluing the opposite faces of a cube by three translations. Then 
the wave numbers of the square-integrable eigenmodes of the Laplacian are restricted to the 
discrete values ki = 2imi/L 1 (i = 1, 2, 3) where n^'s run over all integers. The angular power 
spectrum is written as 



M^{k)Fl 



k,<> 



kl 



1 pio d$ 

-$(v*)ji(HVo-V*))+2 Jdv-^ji(KVo-v)), (5) 



where V$(k) is the initial power spectrum for $ and k = ykf + fc| + k\. From now on we 
assume the scale-invariant Harrison- Zeldovich spectrum (V$(k) = const.) as the initial power 
spectrum. 
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FIG. 1. Suppression in large-angle power for 3-torus models with or without the cosmological 
constant. A large-angle suppression in AT) = + 1)Ci/(2it) (all the plotted values are nor- 
malised by AT20 with infinite volume) occurs for the "standard" 3-torus model (f2 m , J^a) = (1-0, 0) 
at I < l cu t ~ lixR^jL — 1 while such a prominent suppression is not observed for the 3-torus model 
with (S! m , f2A)=(0.1,0.9). The transfer function T/(/c) describes how each /c-mode contributes to 
the power for a particular angular scale I. The corresponding first eigenvalues k\ are plotted as 
horizontal lines(lower figures). The unit of k is Hq. 



Now let us estimate the angular scale below which the power spectrum is suppressed 
owing to the mode-cutoff k cut = 2tt/L. First of all we consider the transfer function Ti(k) 
which is defined by 

that describes how the spatial information is included in the angular power. For a given 
k, the contribution to the angular power on smallest angular scales comes from the OSW 



5 



effect where the transfer function can be written in terms of the spherical Bessel function 
as TP SW [k) = ji(kr) ). Because jz(x)'s have the first peaks at x ~ 1 + 1, the angular cutoff 
l cut is determined by l cut = 2nri /L — 1. For instance, l cut ~ 5 for L = r] Q = 2i/ ~ 1 . Thus 
the large-angle suppression scale is determined by the largest fluctuation scale at the last 
scattering. On smaller angular scales / > l cut , the second peak in the power corresponds 
to the fluctuation scale of the second eigenmode at the last scattering. This behavior is 
analogous to the acoustic oscillation where the oscillation scale is determined by the sound 
horizon at the last scattering. On angular scales larger than l cut one observes another 
oscillation feature in the power for models with a smaller cell where k cut rjo is sufficiently 
large which is apparently determined by the behavior of the first eigenmode. From the 
asymptotic form 

ji(x) ~ -sm(x-l7r/2), x» ^ ^ % (7) 
x 2 

one notices that the angular scale of the oscillation for asymptotic values k cut rjQ » 1(1 + 1) /2 
is Al = 2. For intermediate values of k cu tf]o, Al take much larger values (see figure |l|). 

So far, we have studied the effect of the non-trivial topology on the OSW contribution 
only which is sufficient for constraining the topology of the "standard" 3-torus model with 
Q t ot — 1- However, for low-matter-density models, one cannot ignore the (late) ISW con- 
tribution which is generated by the decay of the gravitational potential at the A-dominant 
epoch. The crucial point is that they are produced at later time 1 + z ~ well 
after the last scattering. Although it is impossible to generate fluctuations beyond the size 
of the cell (in 3-dimensional sense), that does not necessarily mean that any fluctuations on 
large angular scales (in 2-dimensional sense) cannot be produced. Suppose that a fluctuation 
is produced at a point nearer to us, then the corresponding angular scale becomes large if 
the background geometry is flat or hyperbolic. Therefore, one can expect that the suppres- 
sion on large angular scales owing to the mode-cutoff is not stringent for low-matter-density 
models. As shown in figure [I], the angular powers for a model with (Q\,Q m ) = (0.9,0.1) 
are almost flat. In contrast to the "standard" model, the transfer function for low- matter- 
density models distributes in a broad range of k that implies the additional late production 
of the fluctuations which contribute to the angular power. Surprisingly, in the low-matter- 
density models with small volume, the slight excess power due to the ISW effect is cancelled 
out by the moderate suppression owing to the mode-cutoff which leads to a flat spectrum. 
However, as observed in the "standard" 3-torus model, the power spectra have prominent 
oscillating features. The oscillation scale for I < l cut is determined by the first eigenmode. 
The peaks in the angular power correspond to the first SW ridge and the first and the sec- 
ond and other ISW ridges. Is such an oscillating feature already ruled out by the current 
observation? We will see the results of our Bayesian analyses for testing the goodness-of-fit 
to the COBE data in the next section. 

IV. BAYESIAN ANALYSIS 

In general, the covariance in the temperature at pixel i and pixel j in the sky map can 
be written as 
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Mij =< TiTj >=Y,< a imai'm> > WiW v Y lm (n i )Y Vm ,(n j )+ < N,Nj > (8) 
i 

where a\ m is an expansion coefficient with respect to a spherical harmonic Yi m , <> denotes an 
ensemble average taken over all initial conditions, positions and orientations of the observer, 
Tj represents the temperature at pixel i, is the experimental window function that 
includes the effect of beam-smoothing and finite pixel size, lii denotes the unit vector towards 
the center of pixel i and < iVjiVj > represents the noise covariance between pixel i and pixel 
j. If the temperature fluctuations form a homogeneous and isotropic random Gaussian field 
then the covariance matrix can be written in terms of the power spectrum Ci as 

= i- J2(2l + ^WfQP^rii ■ n 3 )+ < > (9) 

where Pi is the Legendre function. Then the probability distribution function of the pixel 
temperature T for the Gaussian field is 



f(T\Ci) = ; ^ exp( -f T -M^iCi) -f |, 

(2tt) w / 2 det 1/2 M{d) V 2 / 



(10) 



where N is the number of pixels. Bayes's theorem states that the probability distribution 
function of a set of parameters A given the data T is 

f(A\T) oc f(f\A)f(A). (11) 

If we assume a uniform prior distribution, i.e. taking f(A) to be constant, the probability 
distribution function of a power spectrum C\ is then 

MCAT) oc exp( -f T -M-Hd) -f ). (12) 

1 ' ^ det 1/2 M(Q) P V 2 ) 



Before applying the method to the flat 3-torus models one must be aware that a flat 
3-torus we are considering is globally anisotropic although it is globally homogeneous. In 
contrast to the standard infinite models, the fluctuations form an anisotropic Gaussian field 
for a fixed orientation if the initial fluctuation is Gaussian. In other words, for a given I, a 
set of aim's {—I <m<l) are not 21 + 1 independent random numbers. In order to see this, 
we write a plane wave in terms of eigenmodes in the spherical coordinates, 

e Jkx = £ b M m ji(kx)Y lm (n), b klm = Y lm (k), (13) 

lm 

where k denotes the unit vector in the direction of k. Then we have 

a im = Y,H k ) himFki (14) 

k 

where is given by (0). Because a>i m is linear in $(&), it is Gaussian. However, they are not 
independent since &fc/ m 's are proportional to spherical harmonics. The statistical isotropy 
is recovered iff the size L of the cube becomes infinite where k takes the whole values for 
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a given k. If one marginalises the likelihood with respect to the 5*0(3) transformation, or 
equivalently, the orientation of the observer, the distribution function of ai m becomes non- 
Gaussian since ai m is written in terms of a sum of a product of a Gaussian variable times 
a variable that is determined by a spherical harmonic 3^ TO (k) where k is a random variable 
which obeys a uniform distribution. Thus in principle, in order to compare a whole set of 
fluctuation patterns over the isotropic ensemble one should use the likelihood function which 
is different from the Gaussian one. 

Nevertheless, we first carry out a Bayesian analysis using the Gaussian likelihood function 
(|P2|) which depends only on C/'s (which is 5*0(3) invariant) in order to estimate the effect of 
the mode-cutoff imposed by the periodic boundary conditions. In the following analysis, we 
use the inverse-noise-variance-weighted average map of the 53A,53B,90A and 90B COBE- 
DMR channels. To remove the emission from the galactic plane, we use the extended 
galactic cut (in galactic coordinates) |p2| . After the galactic cut, best- fit monopole and dipole 
are removed using the least-square method. To achieve efficient analysis in computation, 
we compress the data at "resolution 6" (2.6°) 2 pixels into one at "resolution 5" (5.2°) 2 
pixels for which there are 1536 pixels in the celestial sphere and 924 pixels surviving the 
extended galactic cut. The window function is given by Wi = GiFi where Fi are the Legendre 
coefficients for the DMR beam pattern and Gi are the Legendre coefficients for a circular 
top-hat function with area equal to the pixel area which account for the pixel smoothing 
effect (which is necessary for "resolution 5" pixels since the COBE-DMR beam FWHM 
is comparable to the pixel size) |p4}| . To account for the fact that we do not have useful 



information about monopole and dipole anisotropy, the likelihood must be integrated over 
Cq and C\ in principle. However, in practice we set Co = C\ = lOOmi^ 2 which renders 
the likelihood insensitive to monopole and dipole moments of several mK. We also assume 
that the noise in the pixels is uncorrelated from pixel to pixel which is found to be a good 
approximation |25|. 

From figure |2| one can see that the constraint on the size of the cell is less stringent for 
the low-matter-density model as expected from the shape of the power spectrum. The effect 
of the suppression of the large-angle power is not prominent unless the cell size is sufficiently 
smaller than the observable region (L < OSH' 1 = 0.08 x 2R* for Q m = 0.1). However, 
the likelihood function varies rapidly as the size increases since the power is jagged in I. 
Unfortunately, the peaks in the power at / ~ 6 and / ~ 11 which correspond to the first and 
the second ISW ridge give a bad fit to the COBE data for L ~ 1.7H Q = 0.17 x 2R*. However, 
the parameter range which gives a good fit to the data is wider for the low-matter-density 
model. It should be emphasised that for both flat models there is a parameter region in 
which the fit is much better than the infinite counterpart. For example, the peaks at / ~ 4 
and I ~ 9 in the angular power of the low-matter-density model with L ~ 2.8HQ 1 give a 
much better fit to the COBE data than the infinite counterpart. In fact the quadrapole 
component in the COBE data is very low and the angular power is peaked at / ~ 4. For 
infinite flat A-models with a scale-invariant initial spectrum such features can be a problem 
since the ISW contribution gives an excess power on large angular scales. 
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FIG. 2. The likelihoods of the 3-torus models with (fi m ,fi/y) = (1.0,0) (diamond) and (0.1,0.9) 
(box) relative to the infinite models with the same density parameters are plotted. The likelihoods 
are marginalised over the quadrapole normalisation Q = (5C2/(47r)) 1//2 . An isotropic Gaussian 
approximation has been used for computing the likelihoods. The COBE-DMR data is compressed 
to 924 pixels at "resolution 5". 



Next, we carry out a full Bayesian analysis in which all the elements of < ai m ai' m > > 
are included. Because of the limit in the CPU power we further compress the data at 
"resolution 5" pixels to "resolution 3" (20. 4°) 2 pixels in galactic coordinates for which there 
are 60 pixels surviving the extended galactic cut. Although the information on smaller 
angular scales / > 10 is lost, we expect that they still provide us a sufficient information 
for discriminating the effect of the non-trivial topology since it is manifest on large angular 
scales. The computation has been done in a similar manner as previous analysis except for 
the covariance matrix for which we use (II) instead of (^). The likelihoods A are computed 
for a total of 2000 random orientations for each model using a vector parallel supercomputer 
vpp 800. The approximated likelihoods which depend on only the power spectra are also 
computed for comparison. 
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FIG. 3. The likelihoods of the 3-torus models relative to the infinite models with the same den- 
sity parameters marginalised over 2000 orientations (star) and that using only the power spectrum 
(diamond). The inverse- noise- variance- weighted average map of the COBE-DMR data is com- 
pressed to 60 pixels (resolution 3). All the likelihoods are also marginalised over the quadrapole 
normalisation Q = (5C2/(47r)) 1 / 2 . 



As shown in figure 0, the discrepancy between the likelihood marginalised over the ori- 
entation of the observer and the approximated likelihood is prominent for models with small 
volume in which the effect of the non-trivial topology is significant. Let us estimate the 
size of the cell for which the effect becomes insignificant. For the "standard" model with 
Q m = 1.0, the volume of the observable region is 47ri?^/3~34 which gives the critical scale 
L c ~3.2if " 1 for which the volume of the cell is comparable to the volume of the observable 
region at present. For the low- matter-density model a significant amount of large- angle 
fluctuations are produced at the A-dominant epoch z ~ 1. Therefore one should compare 
the volume of the sphere with radius r\{z = 0) — r](z— 1) ~0.9i/ ~ 1 to the volume of the cell 
instead of the observable region which gives L c ~ 1.5HQ 1 . We can see from figure that 
these estimates well agree with the numerical result. 

We can give two explanations for the discrepancy although they are related each other. 
One is the non-Gaussianity in the fluctuations and another one is the correlations between 
ai m 's owing to the global anisotropy in the background geometry. 

Suppose a random variable Z = XY in which X obeys a distribution function E(X) 
which is even and Y obeys a distribution function F(Y). Then the distribution function G 
of Z is given by 

G(Z) = J ^pF(Y)dY, (15) 

which is apparently even. Because the fluctuations are written in terms of a sum of products 
of Gaussian variable <E>(fc) (with zero average) times a non-Gaussian variable Y/ m (k), the 
skewness in the distribution function of a\ m marginalised over the orientation is zero although 
the kurtosis is non-zero. 

The correlations in a/ m 's are the consequence of the gap between the degree of freedom 
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of (k) and (l,m). The degeneracy number in a fc-mode (=the number of the direction k) is 
much less than the number of relevant "quantum numbers" (/,m)Q if the scale (=27r//c) is 
comparable to the length of the side L. Taking an ensemble average over the initial condition 
we have 

87T 3 

< ai m a Vm > >= T573^(^)^m(k)^' OT '(k)F fc iF W /, (16) 

where Fki is given by ([5|). The sum does not vanish in general when (l,m) ^ (l',m') which 
is the consequence of the global anisotropy in the background geometry. However, if one 
takes an average of ([RJ) over k, one finds that all off-diagonal elements vanish because of 
the orthogonality of the spherical harmonics. Similarly one can consider 4-point correlations 
< CLi imi ai 2m2 ai 3m3 ai 4m4 >. In this case, all the off-diagonal elements (k,mi) ^ (lj,rrij) do not 
necessarily vanish even if one takes an average over k. Thus the non-Gaussianity for the flat 
3-torus models contrasts sharply with the one for the compact hyperbolic models in which 
the pseudo-random Gaussian property of the expansion coefficients bki m (which are obtained 



by expanding the eigenmode in terms of eigenmodes in the universal covering space) |26|j27 
renders off-diagonal elements always vanish if an average is taken over the position of the 
observer (although the kurtosis is non-zero). The difference can be attributed to the property 
of eigenmodes. In a less rigorous manner the property of the eigenmodes which are projected 
onto a sphere with large radius can be stated as follows: Projected eigenmodes of compact 
hyperbolic spaces are "chaotic" whereas those of compact flat spaces are "regular". 

As for the constraint on the size of the cell, we set a slightly severe condition for the 
low-matter-density model since the likelihood of the infinite counterpart is ~ 10 _1 of that 
of the infinite "standard" model due to a slight boost on the large angular scales caused by 
the ISW effect. Together with the previous analysis using the data on the "resolution 5" 
pixels, the conditions on the relative likelihood log 10 (A/A oo ) > —2 and log 10 (A/A oo ) > — 1 
yield L > 1.6# _1 and L > 2.2# ~ 1 for the "standard" 3-torus model (fi m , fi A ) = (1.0, 0) and 
the low-matter-density 3-torus model (ft m , Ra.) = (0.1,0.9), respectively. Here Aoo denotes 
the likelihood of the infinite counterpart with the same density parameters. The maximum 
number N of images of the cell within the observable region at present is 8 and 49 for 
the former and the latter, respectively. Note that the constraint on the "standard" 3-torus 
model is consistent with the previous result f3l.fl . 



2 Here we assumed that for each degenerated mode specified by k, are independent Gaussian 
variables. 
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V. SUMMARY 



In this paper the CMB anisotropy in a flat 3-torus model with or without the cosmo- 
logical constant has been investigated. Using the COBE-DMR data, we have done a full 
Bayesian analysis incorporating the effect of anisotropic correlation. It has turned out that 
the constraint on the low-matter-density model is less stringent compared to the "standard" 
model with Q m = 1.0. The reason is that the large-angle fluctuations (on COBE scales) are 
produced at late time well after the last scattering. The physical size of these fluctuations 
is of the order of the size of the cell or less. Hence the effect of the non-trivial topology 
becomes insignificant. We expect that the result does not significantly change for other 
compact flat models with different topology || or with different shape of the cell since the 
physical effect does not change. 

We have seen that the analysis using the angular power spectrum is not enough since the 
background geometry is globally anisotropic. Even if the power spectrum well agrees with 
the data, it does not guarantee the validity of the model since the power spectrum has infor- 
mation of only isotropic components in the 2-point correlations. If the background geometry 
is globally anisotropic then the fluctuations form an anisotropic Gaussian field for a given 
orientation. On the other hand, the fluctuations become non-Gaussian if one marginalises 
the distribution function over the orientation ||28|| . For locally Friedmann- Robertson- Walker 
models which are spatially compact, we expect to observe zero skewness but non-zero kur- 
tosis provided that the initial fluctuation is Gaussian. 

However, for low-matter density models, anisotropy in the statistically averaged correla- 
tion on large angular scale cannot be so large since the physical size of observed temperature 
fluctuations that has been prodeced well after the last scatteing is much smaller than the 
actual size of the space(i.e. topological identification scale). 

In the case of flat topology, the discrete eigenmodes have "regular" features which lead to 
significant correlations in a/ m 's. Even if marginalised over the orientation, the correlations 
do not completely disappear. This contrasts with the case of hyperbolic topology in which 
the discrete eigenmodes are "chaotic" and the correlations in a/ m 's disappear if one takes an 
average over the position of the observer. 

In order that we would be able to observe the periodic structure in the fluctuations 
for the 3-torus models, we need to have L/(2i?*) < 1 where .R* is the comoving radius of 
the last scattering surface. Because the obtained constraints give L/(2i?*) > 0.40,0.22 for 
(fl m , Q\) = (1.0, 0), (0.1, 0.9), respectively, we still have a great chance of the first detection 
of the non-trivial topology of the universe by the future satellite missions such as MAP and 
Planck which can survey the CMB in the full sky with high resolution by measuring some 
specific signatures (e.g. the circles test El or the non-Gaussianity P7J). 
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